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ABSTRACT 

The first paper in this series explored the effects of altering the chemical mixture of the stellar 
population on an element by element basis on stellar evolutionary tracks and isochrones to the end of 
the red giant branch. This paper extends the discussion by incorporating the fully consistent synthetic 
stellar spectra with those isochrone models in predicting integrated colors, Lick indices, and synthetic 
spectra. Older populations display element ratio effects in their spectra at higher amplitude than 
younger populations. In addition, spectral effects in the photospheres of stars tend to dominate over 
effects from isochrone temperatures and lifetimes, but, further, the isochrone-based effects that are 
present tend to fall along the age-metallicity degeneracy vector, while the direct stellar spectral effects 
usually show considerable orthogonality. 

Subject headings: stars: abundances — stars: atmosphere — stars: evolution — stars: fundamental 
parameters — globular clusters: general — galaxies: abundances — galaxies: stellar 
content 



1. INTRODUCTION 

Little is known about the influence of individual ele- 
ment abundances on the integrated flux of a stellar pop- 
ulation. As present, it is impossible to derive an accurate 
age within 10% from the integrated light of a cluster- like 
single-age and single-abundance stellar population, and 
the primary r eason is the com plication due to abundance 
ratio effects (|Worthevlll998D . The secondary reason is 
that the input ingredients and choices made in even the 
most sophisticated stellar evolution calculations induce 
scatter in the results (Char iot et al.| [r996). The tertiary 
reason is the set of difficulties associated with stellar flux 
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knowledge, such as colors, line strengths, or spectra, that 
are needed at each stellar evolutionary phase to represent 
the component stars. Additional uncertainties exist, such 
as dust extinction, stellar rotation and activity, the blue 
straggler frequency, and the effects of close bi naries. 

Paper I in this series, iDotter et al.l (|2007bf ). explored 
the effects of 12 chemical mixtures on the L, Tcff, and life- 
times of stars along stellar evolutionary isochrones and 
upon the opacities needed to calculate the stellar models. 
The mixtures explored were solar, a-element enhanced, 
and ten cases where only one element at a time was en- 
hanced, for elements C, N, O, Ne, Mg, Si, S, Ca, Ti, 
and Fe. The mixtures were re-scaled so that the mass 
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fraction of heavy elements, Z, was constant. 

The conclusions from that paper could be summarized 
by splitting the elements that were investigated into three 
categories. "Displacers" include O, C, N, and Ne. These 
are elements that are abundant but supply less opac- 
ity per unit mass than heavier elements. At fixed Z, 
boosting a displacer element will therefore decrease the 
opacity. This leads to shorter stellar lifetimes and hotter 
stars. "Boosters" include Mg and Si. These elements are 
good opacity sources, so that boosting a booster will in- 
crease opacity and make cooler stars that live somewhat 
longer. "Oddball" elements defy these trends, and in- 
clude Ca and Ti, both of which make cooler stars, that, 
nevertheless, have shorter lifetimes. Another oddball is 
S, which decreases low-temperature opacity but increases 
high-temperature opacity so that its ultimate effect on 
temperatures and lifetimes is small. An element worth 
special mention is Fe, because its high-temperature opac- 
ity has considerable structure. One of its opacity peaks 
corresponds to the temperatures characteristic of the 
outer edge of the convective core that develops in stars 
slightly more massive than the sun. Increasing the Fe 
abundance therefore strongly couples to the onset of con- 
vective core overshooting effects, and leads to a strong lu- 
minosity effect in the region of the main sequence turnoff. 
The temperature, luminosity, and lifetime effects are il- 
lustrated in Paper I, but further consequences will be 
illuminated in this paper. 

Besides isochrones, the other major ingredient in any 
population synthesis effort is some representation of stel- 
lar flux for each star in the isochrone. The effects of indi- 
vidual element abundances on the spectra of stars have 
been studied for years, and enormous progress has been 
made in finding the structure of stellar atmospheres and 
calculating the emergent flux. The two most-cited pro- 
grams fo r constructing model atmosp heres of stars are 
ATLAS (iKuruczlflQTni i and MARCS (iGiistafsson et all 
fl975h and we use results from these code in this work. 
For sp ectral synthesis of the emerge nt flux we use SYN- 
THE (iKur ucz' 1970^, FAN TOM dC oelho etal] [20051 : 
ICavrel et al. 1991), and SSG (jBell et al. 1994). For ordi- 
nary population synthesis, empirical spectra or colors can 
be used (.Vazdekis.1999) . but for investigating element- 
by-element effects, using synthetic spectra is clearly the 
way forward. Investigations into the effects of individual 
elements on stellar spectra, but with the intent of apply- 
ing th e resu lts to galaxy spectr a, include iTripicco fc Bell 
(| 1995( 1 and iKorn et all (|2005( l. who gauged the effects 
of ten individual elements on th e Lick system of 25 
pseudo-equivalent width in d ices (jWorthe v et al.l 119941 : 
iWorthev fc Ottaviani|[l997l ). iServen et al . '(2005) inves- 
tigated 24 elements in spectra that ran from 3500 A to 
9000 A with velocity smoothing appropriate for dynam- 
ically hot systems such as elliptical galaxies. 

The present paper combines the new isochrones de- 
scribed in Paper I with greatly extended and updated 
synthetic spectra in order to obtain ab initio popula- 
tion synthesis models as a function of individual element 

^ Added to these 12 mixtures were three at non-constant Z and 
variable X and Y but constant [Fe/H] with 0.2 dex more Carbon, 
0.3 de x more Nitrogen and 0.3 dex more Oxygen (see lDotter et aU 
I2007bl 'l. 



abundances^. All parts of the models, including high 
temperature and low temperature opacities, stellar evo- 
lutionary models, and synthetic spectra, include the al- 
tered abundances in the same way. At present the grid 
includes only solar Z , but calculations are underway to 
extend to many different abundances. This paper de- 
scribes the new spectra and their color index results in § 
2, Lick index results in § 3, concluding with a discussion 
and summary section. 

2. DESCRIPTION OF NEW MODEL INGREDIENTS 

The new isochro nes were described in Paper I 
(jPotter et al.l l2007bl ). and we refer the reader to that 
paper for details. Twelve chemical mixtures were ex- 
plored self-consistently with both high-temperature and 
low-temperature opacities adjusted properly, and evolu- 
tion is complete only to the end of the red giant branch. 
It should be emphasized that the models are thus incom- 
plete and should not be used blindly when comparing to 
real stellar populations until the helium-burning phases 
are properly incorporated'^. This is planned for the ongo- 
ing improvement of these models. In careful, differential 
ways, the present models give us important pointers to 
the behavior of stellar populations with variable element 
mixtures, but the absolute values of colors and indices 
are not to be trusted. The mixtures explored were solar, 
a-element enhanced, and ten cases where only one ele- 
ment at a time was enhanced (C, N, O, Ne, Mg, Si, S, 
Ca, Ti, and Fe). The mixtures were re-scaled so that the 
mass fraction of heavy elements, Z, was held constant. 

2.1. Synthetic Spectra 

New for the present work is a collection of synthetic 
spectra from three different sources. Spectra for the 
coolest stars, at 3000 and 3190 K for log g — 0, were 
calculated using the FANTOM synthesis code and iPle j 
([1992 ) model atmospheres. At 3500 and 3750 K for log 
5 = and 5 spectra were calculated using the FANTOM 
synthesis co de and ATLAS model atmospheres exactly as 
described in lCoelho et"all (j2005f ). from 3000 A to 10000 
A in 0.02-A steps. Spectra for the medium-temperature 
stars were calculated using the SSG synthesis code and 
MARCS model atmospheres from 3000 A to 10000 A in 
0.01-A steps. There were 9 temperatures in this range, 
from 4000 K to 6000 K, in 250 K steps. Gravities of log 
5 = 4.5 and 2.5 were used for all nine temperatures, but 
log g = 0.5 was included for the three coolest. Spectra 
for hot stars, at 7000 K, 8000 K, 10000 K, and 20000 K, 
were calculated using the SYNTHE synthesis code start- 
ing from ATLAS model atmospheres in high resolution 
with logarithmic wavelength spacings. Gravities of logg 
= 4.5 and 2.5 were used except the 20000 K models, 
where 4.5 and 3.0 were used. The lines lists were not 
completely homogeneous between the three regimes, al- 
though they share much in common. The hot stars were 
computed with R. Kurucz (kurucz.harvard.edu) line lists, 
the medium-temperature stars included custom modifi- 
cations by R. Bell, M. Tripicco, and M. Houdashelt, and 

^ ICoelho et al.l l|2007l ) did a similar investigation with a flat en- 
hancement of all the a-elements over a range of metallicity. 

^ The contribution from the later stellar evolutionary phases to 
the spectral indices can be inferred from Figure 8 in Coelho et al. 
(2007). 
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the cool stars benefitted from the TiO re-scalings de- 
scribed in ICoelho et al.l ()2005l ). 

These spectra were then rebinned to a common wave- 
length range of 3000 A to 10000 A and a common lin- 
ear wavelength binning of 0.01 A per flux point for fu- 
ture high-resolution studies, and then rebinned again to 
a common linear wavelength binning of 0.5 A per flux 
point for more convenient use for working with colors 
and spectral indices'*. 

For each star in the grid of spectra described above, 
multiple realizations of the same spectrum were cal- 
culated with single-element abundance adjustments, so 
that the effects of C, N, O, Ne, Mg, Si, S, Ca, Ti, and Fe 
could be individually gauged. For neon enhancement, we 
use the scaled-solar spectra because of the small effect on 
the optical stellar spectrum from neon. Any isochrone- 
level temperature, luminosity, or stellar lifetime effects 
that we depicted in Paper I from neon are, however, in- 
cluded. Thus, at solar metallicity, 35 gravity and temper- 
ature combinations x 10 element adjustments (counting 
the non-adjusted solar abundance spectrum) makes a to- 
tal of 350 spectra. 

Another convenience is the definition of i?, which 
stands for "generic heavy element." Usually, one refers 
to relative abundance changes relative to iron, such as 
[0/Fe] or [Mg/Fe]. However, this stops making sense in 
the case of wanting to change Fe with respect to other 
elements; i.e., using [Fe/Fe] leads to nonsense. So we use 
the "generic heavy element" R to stand for "all heavy el- 
ements that remain scaled with the solar mixture, where 
all other element tweaks must be specified." With that 
notation, [Fe/R] is perfectly acceptable, and the specifi- 
cation of constant [R/H] means that all elements except 
the one (or, more generically, ones) under consideration 
are held in solar lockstep. 

2.2. Synthetic Spectra Accuracy 

An interpolation routine was constructed so that spec- 
tra of arbitrary temperature, gravity, and abundance 
mixture could be produced on demand. Linear inter- 
polation in the log of the flux seemed to give the most 
predictive results. Interpolation was done for all input 
variables: temperature, log 5, and the whole set of abun- 
dance parameters. Tests in which one star in the grid 
was compared with its interpolated version calculated 
from the flanking temperature grid points indicate that, 
within the dense 4000-6000 K part of the grid, interpo- 
lations yield better than 1% accuracy at any given flux 
point, since the spectra change more with temperature 
than with the other parameters. Give n the nonli near 
tempe rature response of many lines, e.g..|Gray fc Brown 
^mh . this is encouragingly good. Doubtless, in the fu- 
ture, we will demand a denser temperature grid, but 1% 
accuracy is better than we need for the present study. 

Toward the more fundamental question of how well 
the synthetic spectra match real spectra, there i s littl e 
we can add to disc u ssions in iMartins fc Coelhol (|2007|): 
ICoelho et all (l2005h: iKm-n et al.l (|2005[ ): iTripicco fcEell 
(|1995( 1. IServen et al.l (|2005( ) did discover some too-strong 



lines due to chromium (that is, mistakes in the line list) 
via stellar comparisons, but these lines do not affect the 
present conclusions at all. More worrisome is silicon be- 
cause the features that show Si effects are mostly SiH 
lines in the blue, but these line lists may be immature 
(R. Peterson, private communication). All conclusions 
regarding Si should be regarded with great suspicion un- 
til this question is satisfactorily resolved. 

In Figures 1 and 2, we compare the high resolu- 
tion spectra of the Sun and Arcturus with the grid- 
interpolated synthetic spectra near the regions of the H/? 
and the Mg b lines, respectively. For the solar parameter 
we employ T^ff = 5777 K, log.g = 4.44, and [Fe/H] = 0, 
and for Arcturus Tejf = 4290 K, logg = 1.9, and [Fe/H] 
= -0.7 with [g/Fe] = -1-0.4 dex (e.g.. iPeterson et al.l 
19931: iGriffin fc Lvnas-Gra^ [l999l: iCarretta et al.l l200l 



* We have used synthetic stellar spectra with an elemental en- 
hancement by 0.3 dex in this study. As discussed in Tripicco & 
Bell (1995), however, a 0.15 dex enhancement is used for carbon 
in order to avoid the carbon star abundance regime where carbon 
atoms outnumber oxygen atoms. 



Fulbright et al.l 120071 : iKoch fc McWilHaml 120081) . The 

top panels are the solar spectra and the bottom panels 
are Arcturus. The thicker (black) lines are the observa- 
tions and the thinner (red) ones are the synthetic. The 
degraded low-resolution spectra, Gaussian-smoothed to 
8 A FWHM or about 200 km s"* — somewhat better 
than the best Lick-system resolution, are overlaid on top 
of the high-resolution spectra. We investigate the in- 
dividual element effects over a rather broad wavelength 
region for the Lick indices and broadband colors in this 
paper. In this context, it is useful to find that the no- 
table mismatches between the synthetic and the observed 
spectra in the high-resolution comparisons are smeared 
out and become subtle at low resolution^. 

Table [T] presents some comparisons between synthetic 
colors and empirical color behavior for the stars in the 
temperature range 4000 to 6000 K. Columns 2-5 refer 
to synthetic spectra, and columns 6- 8 refer to the em- 
pirical calibration of IWorthev fc Lei (|2008|) . The syn- 
th etic spec t ra we re converted to UBVRcIc colors as 
in I WorthevI (jl994f ) using Bessell's (1990) filter responses 
zeroed to Vega's colors. There is at least a 10 mmag 
uncertainty just from that procedure, and probably at 
least 20 mmag for U — B. The "dwarf is log g = 4.5, 
and the "giant" is logg = 2.5. Columns 2 and 6 present 
the colors for the star listed, synthetic and empirical, re- 
spectively. The most serious mismatches are the U — B 
color for the 6000 K dwarf and for the 4000 K and 5000 K 
giants. These facts suggest some room for improvement 
in model atmospheres and line lists (especially at shorter 
wavelengths) . 

The remaining columns present color changes induced 
by either an abundance increase of 0.3 dex or a tem- 
perature increase of 250 K (color changes in milli- 
magnitudes). Column 4 is a complete scaled-solar en- 
hancement of 0.3 dex in every element and with new 
atmospheres calculated, while Column 3 is a 0.3-dex en- 
hancement of only the 24 elements that are explicitly 
traced in the spectral library, a superset of those followed 
in the isochronc library.^ Comparison of columns 3 and 
4 indicates that, for the "heart" of the spectral library 

^ Due to assumptions about scattering versus absorption 
in the SSG synthesis code, the bottoms of the saturated 
lines are artificially flattened. Some more comparisons be- 
tween our synthetic spectra and Arcturus can be found at 
http://astro.wsu.edu/hclee/NSSPM_II_Arcturus.html 

° They are C, N, O, Na, Mg, Al, Si, S, CI, K, Ca, Sc, Ti, V, Cr, 
Mn, Fe, Co, Ni, Cu, Zn, Sr, Ba, Eu. 
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between 4000 and 6000 K, the sum of the available one- 
by-one element tweaks approximately equals the scaled- 
solar analog operation. The list of 24 elements in linear 
combination thus appears to approximate the full-blown 
calculation of a spectrum in most cases. This allows the 
approximate calculation of a new spectrum at arbitrary 
composition in an eyeblink rather than many minutes for 
a whole new synthetic spectrum calculation. 

In general, the synthetic colors are within a few hun- 
dredths of a magnitude of the empirical colors and track 
the empirical color responses to abundance and temper- 
ature in an approximate way. The worst outlier is the 
U — B color for the 6000 K giant, in which the color gets 
redder r ather than bluer with in creasing temperature 
(see also iMartins fc Coelhol I2007D . Strong conclusions 
are not possible with this comparison, but the agreement 
is fairly encouraging. 

Tables [2] and [3] show the spectral effects of element- by- 
element enhancements of 0.3 dex (except carbon which 
is increased by 0.15 dex, see footnote #4) on some se- 
lected color indices for the stars in the temperature range 
4000 to 6000 K. In Table [2] the elements were rescaled to 
constant heavy-element fraction Z, while in Table [3] Z 
was allowed to increase. One technical subtlety is that, 
since neon is not tracked in the spectral library, the only 
effects that appear in Tabic [2] for Ne are due to the re- 
scaling, and there are no effects at all in Table [31 Ta- 
ble 3 can be visually explained by taking a look at the 
spectra^. Figure 3 displays one example. Here the flux 
ratio between 0.3 dex Mg-enhanced spectra and that of 
the solar-scaled are shown for a 4000 K giant star. The 
B^V becomes bluer (e.g., more i?-band flux over y-band 
flux) and V — R, V — I become redder (e.g., compara- 
bly more R- and /-band fluxes over F-band flux). The 
U — B also becomes similarly bluer just like B — V be- 
cause the B-band filter response curve goes all the way 
to ~5400 A. Now the effects of Mg that we read from 
Table 3 are clearly understood from this illustration. It 
is confirmed that Mg is the most important a-element 
that influence th e U — B and B — V colors as already 
described in Cassi si et al.l (j2004). 

Comparing corresponding entries in the two tables, es- 
pecially at the 6000 K giant, one sees that oxygen have 
little effect on the spectrum by themselves; it is the ef- 
fect of the decrease of the rest of the elemental abun- 
dances that causes the bulk of the spectral change in the 
constant-Z case. Unsurprisingly, it is also clear from the 
table that cooler stars are more susceptible to element- 
by-element effects than warmer stars, at least at optical 
wavelengths. 

2.3. Integrated-light Color Results 

Table|4]presents color results for the isochrone-summed 
models assumi ng a S alpctcr (1955) initial mass function 
and using the IWorthevI (|l994f l machinery for some dis- 
parate ages of populations. The bolometric correction 
(BC) row includes both isochrone effects and direct stel- 
lar spectral effects with the caveat that our synthetic 
spectra are limited in wavelength coverage, so the fluxes 
at given wavelength s were normalized to the older low- 
resolution IWorthevI ()1994l ) flux library at the same Z. 

They are given at |http://astro.wsu.edu/hclee/| 
NSSPMJLColor.html 



Thus, only effects that affect wavelengths between 3000 
and 10000 A are considered, and there may be additional, 
small effects present that are not accounted for. However, 
in absolute value, all of the BC shifts are quite small, the 
largest being 6% for the case of all a-elements enhanced. 
It is unlikely that element ratio effects will be of major 
concern for the mass estimation of clusters and gal axies 
from their luminosities (see also lCassisi et al.]|2004[ ). 

The other immediate conclusion from Table [4] is 
that older populations show larger spectral effects due 
to element-by-element abundance changes. This is a 
straightforward consequence of the earlier conclusion 
that cooler stars are more sensitive; the light of older 
populations is dominated by stars that are cooler than 
those present in younger ones*. Oxygen and neon tend 
to make the colors bluer, while species that contribute 
more to the lines in the spectrum generally make them 
redder. Exceptions can be easily explained if one knows 
where the lines contribute the most. For instance, mag- 
nesium has most of its absorption near 5100 A, I.e., m 
the V band, so adding Mg makes B — V bluer while it 
makes V — R and V — I redder. The former, the bluer 
B — V \s, also partly due to the increase in the flux in the 
region around 4000 A from the Mg -enhanced spectra a s 
shown in Figure 3 and described in lCassisi et al.l (j2004f ). 
The latter, the redder V — R and V — 1 also reflect the 
cooler red giant branch as illustrated in figure 8 of Paper 
I. Figure 4 displays this Mg-enhanced integrated case at 
5 Gyr (thick line) over the single star case of a 4000 K 
giant (thin line). The populations respond rather like the 
stars do, as can be seen from Figure 4. Based on this, 
one would suspect that isochrone-caused effects are rel- 
atively minor though non-negligible, and this conclusion 
will be confirmed and amplified in the following section. 

3. RESULTS IN LICK INDEX DIAGRAMS 

The effects of element by element enhancement on 
the Lick indices (jWorthev et al.lfl994f ) are described in 
this section^. The synthetic indices are not very ac- 
curate in absolute predictions (c.f. iKorn et all l2005t 
iServen et all [2001 so we employ a differenti al ap- 

J roach in w hich the fitting fu nctions of Worthc v et al.l 
1991 and iWorthev fc Ottavianil (|1997l) are used as 
the zero point, and delta-index information as a func- 
tion of element ratio is incorporated via measuring the 
synthetic spectral library. This procedu re is similar 
to tha t of previous inve stigations (e.g.. iTrager et"aLl 
2000a"l?; Proctor & Sansom 2002; T homas et all l2003t 
Lee h Worthev. ,2005 : Schiavon 2007,) but more sophis- 
ticated since an entire grid of delta-index information 
was used. That is, 350 spectra at solar Z, plus similar 
grids for 4 other Z values for this work, as opposed to 2 
or 3 synthetic stars at solar abundance only for previous 
works. 

To maintain an exact correspondence with Paper I, all 
the elements are 0.3 dex enhanced, except carbon, which 

^ This last statement is vaUd for ages > 1 Gyr. For younger ages, 
IR fluxes from AGB and TP-AGB stars complicate this status (see 
also ILee et al.li2007ai ') . 

® In this study, we mostly describe the elemental effects on the 
Lick indices from the integrated spectra. Prompted by the referee's 
suggestion, however, we have looked into those elemental effects on 
th e Lick indices at the stellar level. Some examples can be found 
at |http://astro.wsu.edu/hclee/NSSPM_II_Lick.html| 
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is enhanced 0.2 dex . In the cases of carbon-, nitrogen- 
, and oxygen-enhancement, we investigate both at fixed 
total metalhcity and at fixed [Fe/H] (and [R/H]). The 
cases at fixed [Fe/H] are denoted with plus sign in the 
figures (e.g., C-I-, N-I-, 0+). Also in the figures, solar- 
scaled solar metallicity predictions are connected by a 
solid line from 1 Gyr to 12 Gyr and the element-enhanced 
cases are marked at 1, 2, 5, and 12 Gyr. 

We have selected Fe5406 as a reference index in most 
plots. Among eight Lick iron indices (Fe4383, Fe4531, 
Fe5015, Fe5270, Fe5335, Fe5406, Fe5709, Fe5782) we pre- 
dict that Fe5406 is insensitive to every element except 
iron (see Figures 10—11 and Table 5 and online spec- 
tra described in footnote #9), making it a convenient 
independent variable. 

As has alread y been mentioned in the literature 
(|Lee et al.ll2007b[ ). for most of cases the isochrone efi^ects 
are relatively minor compare d to the stellar spectral ef- 
fects (see also ISchiavonll2007i ) . For some cases, however, 
isochrone effects are non-negligible and quite i mportant 
(cf. H ff in Figure 11; see also figure 17 in iCoelho et al.l 
[200l . 

It may be worthwhile to mention for clarity that 
"isochrone effects" indicate the temperature, luminosity, 
and stellar lifetime effects from element ratio changes, 
and, in the case of the iron-enhanced mixture, the al- 
tered [Fe/H] value that goes into the empirical fitting 
functions in our experiment. When calculating observ- 
ables, we would still use a scaled-solar-ratio spectral li- 
brary. These isochrones were presented in Paper I. In 
the present Paper II we add the detailed spectral effects 
due to element by element enhancement at the stellar at- 
mosphere/stellar flux level. The "direct stellar spectral 
effects" are the ones that come purely from the emergent 
spectra with isochrones held fixed. 

[Carbon and Nitrogen: CN2]: Lick indices CNi and 
CN2 have identical central bandpasses, but CN2 has a 
narrower blue continuum which makes the CN2 some- 
what less prone to abundance variations other than those 
due to elements C and N. We compare CN2 with G4300 
as carbon, nitrogen, oxygen, silicon and iron are en- 
hanced in Figure [5l The upper panels show that CN2 
is both carbon- and nitrogen-sensitive. Nitrogen's effect 
is more prominent (partly because carbon is enhanced 
only by 0.2 dex compared to the 0.3 dex nitrogen en- 
hancement). The bottom right panel also suggests that 
silicon and iron affect CN2 in a non-negligible way when 
stellar populations become older than 5 Gyr. We can 
tell from Table 5 that the silicon effect is mostly from 
the stellar spectra, while the iron effect is largely from 
the isochrones. The oxygen-enhanced case is displayed 
in the bottom left panel. It shows that simply adding O 
alters CN2 and G4300 by a small amount (increases at 
young ages and decreases at old ages; the latter via con- 
sumption of more C into the CO molecule via molecular 

As footnote #4 says, the carbon-enhanced spectra are gener- 
ated with 0.15 dex carbon enhancement. But the carbon-enhanced 
isochrones that we presented in Paper I are of 0.2 dex carbon 
enhancement. In order to be consistent, for the stellar popula- 
tion synthesis calculations, we extrapolated those 0.15 dex carbon- 
enhanced spectra to 0.2 dex enhancement in order to match the 
carbon-enhanced isochrones. The figure 1 of Paper I, in fact, needs 
to be corrected. The filled box for the case of carbon should be 
located near —0.04 instead of near in terms of [Fe/H]. 



equilibrium balance^"'^), but if Z is held constant both in- 
dices decrease much more due to the displacement of C 
and N to lower abundance because of the 0-enhancement 
and fixed sum. G4300 is N-insensitive. Also, we note 
that CN2 has little sensitivity to Mg-, S-, Ca-, or Ti- 
enhancement. 

[Carbon: G4300 vs. C24668]: In Figure d G4300 
and C24668 are plotted as a function of Fe5406. The 
left panels show C- and Fe-enhanced models, while the 
right panels show the effects of Si and Ti enhancement. 
G4300 and C24668 are known to be good carbon in- 
dicators among Lick indices along with CNi and CN2. 
The left panels show that this is indeed the case. In the 
bottom left panel, however, it can be seen that C24668 
is also highly iron sensitive, as indicated by its former 
name, Fe4668 (Worthey 1994)^^ Contrary to C24668, 
G4300 shows a negligible iron sensitivity. Furthermore, 
the right panels show that C24668 is influenced by Si and 
Ti, whereas G4300 shows only slight Ti sensitivity. Both 
G4300 and C24668 show little sensitivity by N-, Mg-, S-, 
or Ca-enhancement. 

The prospects for disentangling C, N, and O are now 
fairly good, because it appears that there is a lot of sen- 
sitivity among various indices. However, at the moment 
it appears that, with the indices available there will be 
considerable degeneracy among the three quantities. In 
the future, adding NH, CH, and CO features may help 
considerablv fe.g..lYong et al.l l2008t iMartefl et"all l200a 
iMarmol-Oueralto et al.ll2008l ). 

[Calcium: Ca4227 vs. Ca4455]: Figure [7] compares 
Ca4227 and Ca4455 with Fe5406. Carbon-, nitrogen-, 
and oxygen-enhanced cases at fixed [Fe/H] are displayed 
at the left panels, while right panels depict calcium- and 
iron-enhanced cases at fixed Z. It is clear from up- 
per right panel of Figure 7 that Ca4227 is significantly 
boosted with calcium enhancement, the effect increasing 
with age. To a lesser degree, Ca4227 is also affected by 
C, N, O, and Fe. Carbon- and nitrogen-enhancement 
make Ca4227 weaker by 0.3 A, whereas the oxygen- 
enhancement does the opposite mostly due to their ef- 
fects at the blue-continuum (see also Proch aska et al. I 
[2OO5I) . 

Contrary to Ca4227, Ca4455 is hardly influenced by 
any of those elements. According to this study, Ca4455 
is found to be the most element-enhancement-free Lick 
ind ex. This is co nsistent with the previo u s find ings 
by iTripicco fc Bell (|1995l) and iKorn et all (120051 ) al- 
though our presentation is based on a large spectral grid 
weighted by isochrones rather than 3 stars. One concern, 
however, is that two recent data sets of Milky Way globu- 

It is worth reiterating what Schiavon (2007) noted in his 
footnote #5. Carbon-enhancement causes the opposite effect from 
oxygen-enhancement and vice versa. This is because the high dis- 
sociation potential of CO molecule. Therefore, at cooler tempera- 
tures, more carbon translates to more CO, resulting in less oxygen 
and vice versa. 

1^ We find from Table 5 that the Fe-sensitivity of C24668 (also 
CN2 and Fe5015) is mostly an isochrone effect. However, it is, in 
fact, not because of the temperature and/or luminosity changes, 
but because of the [Fe/H] changes that go into the fitting function 
that we use for the index calculation. As one can find from Figure 1 
of Paper I, [Fe/H] = 0.268 for the Fe-enhanced isochrones compared 
to —0.225 for the a-enhanced ones at constant solar metallicity, Z. 
Clearly, future fitting-function work will need to be more meticu- 
lously defined in terms of abundance parameters, and not locked 
to [Fe/H] necessarily. 
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lar clusters, by ICohen. Blakeslee. fc Rvzhov! (|1998f ) and 
iPuzia et al.1 (|2002f ). both significantly disagree with the- 
oretical model predictions and also each other (Figure 3 
of iLee fc Worthcv 2005). According to Tables 1 - 3 of 
iTripicco fc Bell (j 19931 . Ca4455 has the strongest depen- 
dence on the bandpass placement (wavelength shift error) 
among Lick indices. This is because the blue-continuum 
and the index bandpass of Ca4455 share the strong Ca 
absorption line feature near 4455 A and consequently 
cancel out its effect, but only if the wavelength match 
is perfect. This phenomenon at the stellar level with 10 
km/sec velocity dispersion is displayed in Figure 8. Both 
Ca4227 and Ca4455 are rather insensitive to Mg-, Si-, S-, 
or Ti-enhancement. 

[Titanium: Fe4531 and Fe5015]: Fe4531 and Fe5015 
are di splayed in Figure 9 with Fe5406. TB95, Trager 
(|1997t ). KMT05, LW05, and lServen eTall (|2005D indicate 
that these two indices are titanium-sensitive, and we con- 
firm that their titanium sensitivity is indeed strong and 
almost comparable to the iron sensitivity. Fe5406, on the 
other hand, demonstrates little titanium sensitivity. 

[Balmer lines: H/3, Hja, ^Jf, H(5^, HcJi?]: Balmer 
lines are widely used as an age indicator because of 
their nonlinear temperature sensitivity in stars, tracing 
better than many indices the temperature of the main- 
sequence tur noff. However, Xee_fe _ ^yorthev (_20_05.) and 
earher work (IWorthev et al.l 119941 : iThomas et aH l2004 
ICoelho et al.l 120071 ) found that they are also abundance 
sensitive to some degree. In Figures fTOHTSl we look into 
their element by element sensitivity in detail. 

Effects of the individual 10 chemical elements' en- 
hancement on the H/3 and Fe5406 are shown in Figure 
10. The effects of carbon, nitrogen, and oxygen enhance- 
ment are of great importance but they are relatively dif- 
ficult to understand here because our experimental setup 
preserves the total metallicity. In other words, it is not 
straightforward to determine whether we are seeing the 
effects of C-, N-, 0-enhancement or whether we are see- 
ing other elements (such as Mg and Fe) counter-effects 
due to their depression. Hence, in Figure 11, we dis- 
play C-, N-, 0-enhancement cases both at the fixed total 
metallicity and at the fixed [Fe/H] (the plus signs). It is 
seen from left panels of Figure 11 that unlike nitrogen, 
oxygen at young ages and carbon at old ages influence 
H/3. The top and bottom right panels of Figure 11 il- 
lustrates that it is mostly the isochrone-level effects of 
carbon and oxygen enhancement that affect H/3. Fig- 
ure 10 further shows that H/3 is similarly altered by the 
Mg enhancement and in this case it is mostly due to the 
synthetic spectra (see also Table 5 and Figure 12). 

In the lower right panel of Figure 10 (also in the bottom 
panels of Figure 11), a grid with a range of metallicity 
([Fe/H] = -2.0, -1.0, -0.5, 0.0, and 0.5) is displayed at 
1, 2, 5, and 12 Gyr. These add itional calculations a re of 
solar-scaled chemical mixtures (jPotter et al.ll20d7a[ ) . We 
have also depicted a-clements enhanced cases both at the 
fixed Z and at the fixed [Fe/H] (Alpha-f)^^. It is intrigu- 
ing to note that the effect due to the Fc-enhancement by 
0.3 dex is seen over the [Fe/H] = 0.5 grid line. This is be- 

The o-elements enhanced case at the fixed [Fe/H] (Alpha-|- ) 
is [a/Fe] = +0.4 dex at [Fe/H] = using the lDotter et al.l l|2007al V 
The [Fe /H] of the a-elements enhanced case at the fixed Z (Alpha) 
is —0.225 using the Paper I. 



cause these grids only reflect the mere isochrone effects of 
iron variation with solar-scaled chemical mixtures, while 
the 0.3 dex Fe-enhancement case shows both isochrone 
and spectral effects combined. 

We find that the H/3 becomes weaker in the a-elements 
enhanced case at fixed [Fe/H] mostly because of the ef- 
fects of Mg (both isochrone and spectral as one could 
see from Table 5), while it stays nearly unchanged in the 
a-elements enhanced case at the fixed Z mostly because 
of the reflection of depression of Fe. At 12 Gyr, a com- 
parison of integrated spectra uncovers that adding Mg 
brings down the blue continuum, making H/3 weaker 
but adding Fe brings down the central bandpass as well 
as the red continuum, the net effects of which tend to 
cancel out effects from Fe. Alpha-enhancement at fixed 
Z mostly reflects the decrease of Fe due to dilution rather 
than overt, direct spectral effects from a-elements. They 
are displayed in Figures 12 to 14. The solid lines are the 
index bandpass and the dashed lines are the blue- and 
the red-continuum edges. 

Figures 15 and 16, on the contrary, display that both 
H7 and H(5 become mildly ('^ 0.7 A) stronger with Mg- 
enhancement. Furthermore, top panels of Figures 15 and 
16 illustrate that the broader A indices (H7A, H(5a) show 
sensitivity to Fe-enhancement and become significantly 
weaker with increasing age, by up to 3 A at 12 Gyr for 
the case of the H(5a. This is mostly due to spectral ef- 
fects, as can be seen from Table 5. They are illustrated 
in Figures 17 (H7) to 18 (H(5). From this experiment, it 
seems that for the high order Balmer lines, the F indices 
(H7i?, H(5i?) are less prone to abundance changes because 
of their narrower index definition and therefore may pos- 
sibly serve as more robust age indicators. Furthermore, it 
is found that 115 p among Balmer lines has the least sen- 
sitivi ty to C, N, and O (see also ISchiavon et al. 1 12002| . 
Proch aska et al. Il2007| ). 

In Figure 19, we show some selected diagrams for il- 
lustrating element enhancement. In the upper left panel, 
CN2 and G4300 are compared and the effects of carbon 
and nitrogen enhancement are depicted. The plot shows 
that CN2 is sensitive to both carbon and nitrogen. CN2 
index values go up by 0.04 mag at 12 Gyr with 0.3 dex 
nitrogen enhancement here. G4300, on the contrary, is 
primarily a carbon sensitive index. Ca4227 and G4300 
are contrasted in the lower left panel of Figure 19 and the 
effects of sulfur, calcium, and titanium are illustrated. It 
is seen that Ca4227 is significantly affected by calcium 
enhancement whereas G4300 is not much altered by any- 
thing but carbon. In the upper right panel, Mg b and 
Fe5406 are compared and the effects of magnesium, and 
silicon are shown. It is clear that Mg b is notably af- 
fected by the magnesium enhancement. H/? and Fe5406 
are contrasted in the lower right panel and, as we have 
seen in Figures 10 and 13, H/3 is not Fe-enhancement 
sensitive, while Fe5406 is significantly affected by Fe. 

Figure 20 is basically same as Figure 19, but here the 
isochrone effects are decoupled from those of synthetic 
spectra. The lines show the combination of isochrones 
and stellar spectral effects, while the points depict the 
isochrone effects alone at 1, 2, 5, and 12 Gyr. It is seen 
again from Figure[20]that the synthetic spectra make the 
spectral lines stronger and/or weaker, while isochrones 

This was hinted in ITripicco fc Belli | |T995|) . 
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play a comparatively minor role. The displacement of 
points in the bottom right panel is mostly due to the 
different [Fe/H] values that go into the index value cal- 
culations rather than the isochrone changes (see Figure 
1 of Paper I). 

4. SUMMARY AND DISCUSSION 

We have commenced this project to make stellar pop- 
ulation models which incorporate flexible chemistry, so 
that almost any interesting chemical mixture can be in- 
terpolated. Paper I dealt with the effects on the stel- 
lar evolution models, examining the temperature and lu- 
minosity changes due to the altered opacites when ten 
chemical elements are individually tweaked to the end of 
the red giant branch. In this paper we combine those 
isochrone effects with the stellar spectral effects in or- 
der to investigate their mixed effects on the integrated 
spectrophotometric indices as well as in their integrated 
spectra themselves. We again emphasize here that the 
models in this study are incomplete in terms of inclu- 
sion of all stellar evolutionary phases and should not be 
used blindly when comparing to real stellar populations 
until the helium-burning phases are properly incorpo- 
rated. A version with a full range of metallicity and with 
horizontal-branch and asymptotic giant branch stars in- 
cluded is planned. Comparison of our models with ob- 
servations of Virgo cluster galaxies will be presented as 
well (Serven et al. in preparation). 

Within our spectral coverage (3000 A to 10000 A), we 
have investigated the broadband color behaviors in the 
UBVRcIc filter set. As one would expect, older pop- 
ulations show larger spectral effects due to element-by- 
element abundance changes. This mostly reflects tem- 
perature effects. We have also confirmed that Mg is the 
most important a-element that sha pes the U — B an d 
B — V colors as already depicted in Cassi si et al.l (|2004l ). 

From our investigation of Lick indices using the inte- 
grated spectra, we find that (1) CN2 is a useful nitrogen 
indicator once we have good carbon abundances from 
G4300 and C24668, but good silicon (and also titanium 
for the C24668) abundances are also needed, (2) Ca4227 
is a robust calcium indicator with some good constraints 
of C, N, and O, (3) Fe4531 and Fe5015 are very useful ti- 
tanium indicators where an independent iron abundances 
is provided, (4) Mg b and Fe5406 are good magnesium 
and iron indicator, respectively, and (5) the variation of 
individual elements affects the Balmer lines. We defer 
the investigation of NaD and TiOi and Ti02 indices un- 
til we have the Na-enhanced isochrones and can model 
TiO molecular effects with confidence. Below, we illumi- 
nate some points that we have described here with the 
help of full spectra. 

Fe4531 and Fe5015: From Figure 9 we see that both 
Fe4531 and Fe5015 prove to be good titanium indica- 
tors. They are almost equally sensitive to t itanium and 
iron. Figures 21-23 show the SAURON (jBacon et all 
l2001t ) spectral range, 4810 A to 5350 A, that includes 
H/3, Fe5015, Mgi, Mg &, Fe5270 and a part of Fe5335. 
Here the integrated spectra of solar-scaled, Mg-, Fe-, and 
Ti-enhanced cases are shown at 2 Gyr with 300 km/sec 
velocity dispersion normalized at 4750 A. H/3 is not in- 
fluenced much by the enhancement of these elements at 
this age, as Figure 10 also shows. Mg b is mostly sen- 



sitive to Mg, as Figure 21 shows, while Fe5015, Fe5270, 
and Fe5335 are notably sensitive to Fe-enhancement, al- 
though Fe5015 is equally sensitive to titanium, as Figures 
22 and 23 display. 

Moreover, we have found that the Fe-enhanced and 
Ti-enhanced spectra show their centroids at different 
wavelengths in the case of Fe4531, opposite sides from 
the scaled-solar spectra, because of different locations of 
iron and titanium lines in their index bandpass. Fig- 
ure [23] displays this case^^. We have preliminary hints 
from high-S/N, high-resolution spectra of Virgo cluster 
galaxies taken at the Kitt Peak 4-m telescope (Serven et 
al. in preparation) that galaxies with the higher veloc- 
ity dispersion tend to have their centroids at the longer 
wavelengths, indicating higher Ti-enhancements. If con- 
firmed, this would indicate a Ti-a r elation, similar to th e 
well-known Mg-cr relation (see also lMilone et aini2000D . 

Balmer lines: 

H/3: The bottom right panel of Figure 10 shows that, 
in an a-enhanced mixture at the fixed [Fe/H] (Alpha+), 
H/3 becomes weaker. This sounds a bit different from the 
conclusion of LW05 where H/3 is a-element insensitive. 
As we can see from Table 5, the effects that we see from 
the bottom right panel of Figure 10 are mostly isochrone 
effects, especially at young ages. The H/3 index becomes 
weaker because of the decrease of temperature of stellar 
isoc hrones from the g-enhanced case a t fixed [Fe/H] (see 
also lCoelho et al.ll2007l: [Schiavonll2007[ ) . LW05 used only 
solar-scaled isochrones at solar and super-solar metallic- 
ities, although these authors employed modified stellar 
spectra because of the a-enhancement. 

Moreover, it is worth emphasizing that a-enhanced 
isochrones at fixed Z are now no longer significantly 
hotter than the solar-scaled ones. This has been re- 
vealed by Paper I, in its figure 11, after incorporating 
the latest deve l opmen t of low-temperature opacities by 
iFer guson et al.l ()2005[ ). If this behavior is indeed true, 
it wil l certainly affect the conclusions of earlier works 
(e.g .. iThomas fc MarastoiJ |2003[ ) that are based upon 
the 'Salasnich et al. isochrones, which were later 

shown by Weis s et al. I ()2006) to present problems. 

H7 and H(5: Fe-enhancement makes these indices sig- 
nificantly weaker, especially for the broader indices, H'Ja 
and RSa (see also iGraves fc Schiavonl[2008h . Fi gures 17 
and 18 show that this is because of Fe-lines that inhabit 
the continuum region. It is suggested therefore that the 
narrower indices (ajF and HSp) may serve as better age 
indicators in case the stellar systems are of wildly differ- 
ent chemical mixtures. 

Our goals are to narrow down the uncertainties of 
mean age estimations to 10% when derived from a sin- 
gle integrated light spectrum. We believe that we can 
reach these goals as we are able to take into account 
the effects from individual chemical elements. The ul- 
timate goal then obviously would be the understand- 
ing of the chemical enrichment and star formation his- 
tories of galaxies, stepping away from the study of 
mean ages and mean metallicities to begin to tackle 
the more realistic problem of multi-age, multi-metallicity 

This phenomenon at the stellar 
level with higher resolution can be seen at 
http://astro.wsu.odu/hclee/Fe45_4000mp00gp25_ti.pdf and 
http: / /astro. wsu.edu/hclee/i<'e45_4000mp00gp25_te.pdt^ 
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stellar populations as they exist in galaxi es. We believe 
that several effo rts such as SPECKM AP (jOcvirk et al. I 
[20061) ■ MOPED (iPa nter et al. 1120071) . and STARLIGHT 
()Cid Fernandes et a l. 2008) will be even more useful once 
they use sophisticated stellar population models with 
flexible chemistry as inputs. 

We thank Jim Rose, Jon Fulbright, and Andreas Koch 
for useful discussions on the Mg-enhancement of the 



metal-rich cool giants. We are also grateful to the referee, 
Ricardo Schiavon, for a thorough reading which greatly 
helped us to improve the paper. Support for this work 
was provided by the NSF through grant AST-0307487, 
the New Standard Stellar Population Models (NSSPM) 
project. H.-c. Lee also acknowledges JP Blakeslee's sup- 
port via NASA grant. PC acknowledges the support by 
the European Community under a Marie Curie Interna- 
tional Incoming Fellowship. 
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TABLE 1 

Gross Color Behavior With Metallicity Increase 



Color 


Synthetic 


A Color, 


A Color, 


A Color 


Empirical 


A Color, 


A Color, 




Color 


24 Elements- 


All Elements- 


Ted 


Color 


All Elements- 








enhanced 


enhanced 






enhanced 






(mag) 


(mmag) 


(mmag) 


(mmag) 


(mag) 


(mmag) 


(mmag) 



Dwarf, Tgff : 

U- B 
B-V 
V-R 
V -I 

Dwarf, Teff : 

U-B 
B-V 

V-R 
V -I 

Dwarf, Tefi 

U-B 
B-V 
V-R 
V -I 

Giant, T^g - 
U-B 
B-V 
V-R 
V -I 

Giant, T^g - 
U-B 
B-V 
V-R 
V -I 



Giant, 
U-B 



'eff ■ 



B - 

V - 

V - 



V 
R 
I 



-- 4000 K 

1.26 
1.43 
0.89 
1.60 

= 5000 K 

0.61 
0.91 
0.52 
0.98 

= 6000 K 
0.02 
0.58 
0.33 
0.66 

: 4000 K 
1.34 
1.45 
0.83 
1.51 

5000 K 
0.57 
0.93 
0.48 
0.92 

: 6000 K 
0.16 
0.55 
0.31 
0.61 



21 
43 
16 
13 



123 
55 
21 
15 



113 
39 
15 
16 



115 
39 
5 
-4 



159 
52 
18 
9 



79 
48 
12 
12 



81 
35 
24 
24 



111 
45 
19 
14 



95 
37 
11 
5 



137 
57 
21 
20 



147 
60 
12 
-1 



83 
54 
9 
4 



-106 
-137 
-94 
-176 



-179 
-90 
-62 

-103 



-33 
-50 
-28 
-55 



-182 
-166 
-120 
-198 



-159 
-107 
-46 
-80 



19 
-64 
-35 
-69 



1.36 ±0.07 
1.36 ±0.02 
0.82 ±0.01 
1.58 ± 0.02 



0.71 ± 0.07 
0.96 ± 0.02 
0.52 ±0.01 
0.96 ± 0.02 



0.16 ±0.07 
0.61 ±0.02 
0.34 ±0.01 
0.65 ±0.01 



1.69 ±0.07 
1.44 ±0.02 
0.79 ± 0.01 
1.51 ±0.02 



0.77 ±0.07 
0.97 ±0.02 
0.49 ± 0.01 
0.94 ± 0.01 



0.11 ±0.07 
0.55 ± 0.02 
0.31 ±0.01 
0.64 ± 0.01 



58 
22 
3 
-5 



160 
62 
11 
13 



114 
42 
4 
1 



52 
22 
5 
4 



161 
63 
10 
11 



109 
41 
4 
1 



-85 
-81 
-74 
-176 



-169 
-101 
-55 
-97 



-75 
-70 
-33 
-60 



-194 
-103 
-82 
-182 



-194 
-116 
-52 
-87 



-81 
-85 
-36 
-70 



Note. — 1. The ""Dwarf'" is log (; = 4.5. and the "Giant" is logp — 2.5. 2. Columns 2 presents the calculations at [R/H] — 0, where R is the 
generic heavy element (see section 2.2). 3. Column 3 presents the case where 24 elements (C, N, O, Na, Mg, Al, Si, S, 01, K, Ca, Sc, Ti, V, Or, 
Mn, Fe, Co, Ni, Cu, Zn, Sr, Ba, Eu) are enhanced by 0.3 dex. 4. Columns 4 and 7 present the calculations at [R/H] = 0.3 dex. 5. Columns 5 and 
8 present the color changes when the temperature is increased by 250 K, synthetic and empirical, respectively. 
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TABLE 2 

Color Behavior With Enhancements At Constant Z 
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= 5000 K 
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Note. — 



1. All elements scaled individually by -|-0.3 dex, except C, which is increased by -|-0.15 dex. 2. Numbers are in milli-magnitude. 
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TABLE 3 

Color Behavior With Enhancements At Constant [R/H] 
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Note. — 



1. All elements scaled individually by +0.3 dex, except C, which is increased by +0.15 dex. 2. Numbers are in milli-magnitude. 
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TABLE 4 

Model Population Color Behavior With Enhancements At Constant Z 



Synthetic AI AI AI AI AI AI AI AI AI AI AI 

Color Color (C) (N) (O) (Nc) (Mg) (Si) (S) (Ca) (Ti) (Fc) (a) 

(mag) (mrnag) (mmag) (rnmag) (mmag) (mmag) (mmag) (mmag) (mmag) (mmag) (mmag) (inmag) 
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78.6 



Note. — 



1. All elements scaled individually by +0.3 dex, except C, which is increased by +0.2 dex (see footnote ^10). 



TABLE 5 

Model Population Index Behavior With Element Enhancements 
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TABLE 5 — Continued 



Index 


Index 


AI 


AI 


AI 


AI 


AI 


AI 


AI 


AI 


AI 


AI 


AI 


AI 


AI 


AI 


Name 


Value 


(C) 


(C+) 


(N) 


(N+) 


(O) 


(0+) 


(Mg) 


(Si) 


(S) 


(Ca) 


(Ti) 


(Fe) 


(a) 


("+) 



H7F -1.8520 -0.3450 -0.5269 0.0450 -0.0228 0.7050 0.1915 0.1800 0.0320 0.0250 -0.0530 -0.0540 -0.6990 0.9790 0.5160 



Note. — 1. The units of CNi, CN2, Mgi, Mg2, TiOi, Ti02 are mag. That of the rest are A. 2. Index Value in the second column is the Index Value of solar-scaled case. 3. A I = Index Value 
of each element-enhanced case — Index Value of solar-scaled. 4. C+, N+, 0+, and a+ are those element-enhanced cases at fixed [Fe/H], whereas the rest are element-enhanced cases at fixed Z. 
5. C and C+ are 0.2 dex enhanced cases and a+ is 0.4 dex enhanced case, while the rest of cases are 0.3 dex enhanced cases. 6. a-elements are O, Ne, Mg, Si, S, Ca, and Ti in this study. 
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4856 4858 4860 4862 4864 4866 
Wavelength (A) 

Fig. 1. — The synthetic and observed spectra are compared near the region of the H/3 Unes. The top panel is the solar spectra and 
the bottom panel is Arcturus one. The thicker (black) lines are the observations. The degraded low-resolution spectra at 200 km/sec are 
overlaid on top of the high-resolution spectra. It is seen that the notable mismatches between the synthetic and the observed spectra at 
the high-resolution comparisons are smeared out at the low-resolution. 




Fig. 2. — Similar to Fij 



igure 1, but near the region of the Mg b lines. 
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Wavelength (A) 

Fig. 3. — At 4000 K and log g = 2.5, solar-scaled and Mg-enhanced spectra are divided. U BV RI-ha.nd filter response curves are overlaid. 
The blue colors o{ U — B and B — V and red colors oiV — R and V — I from the Mg-cnhanced spectra that we read in Table 3 are visually 
elucidated here. 
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Fig. 4. — Similar to Figure 3, but here the thicker line is the flux ratio between the solar-scaled and the Mg-enhanced spectra of the 5 
Gyr integrated spectra at solar metallicity. It is normalized at 6000 A. The color changes because of the Mg-enhancement that wc read 
from Table 4 can be vividly seen here (see text). 
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Fig. 5. — CN2 is plotted as a function of G4300 from 1 Gyr to 12 Gyr for simple stellar populations at solar metallicity. Carbon-enhanced 
(upper left), nitrogen-enhanced (upper right), oxygen-enhanced (lower left), and silicon-, iron-enhanced cases (lower right) at fixed Z are 
depicted at 1, 2, 5, and 12 Gyr. C-I-, N-I-, and O-f indicate the carbon-, nitrogen-, and oxygen-enhanced cases at fixed [Fe/H], respectively. 
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Fig. 6. — Fe5406 vs. G4300 (top) and C24668 (bottom). Left panels show carbon- and iron-enhanced cases, while right panels display 
silicon- and titanium-enhanced cases. Note that C24668 is significantly affected both by carbon and iron. C24668 is further sensitive to 
silicon- and titanium-enhancement. 
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Fig. 7. — Fe5406 vs. Ca4227 (top) and Ca4455 (bottom). Left panels show carbon-, nitrogen-, and oxygen-enhanced cases at fixed [Fe/H], 
while right panels depict calcium-, iron-enhanced cases at fixed Z. Ca4227 is predominantly calcium-sensitive as well as a contribution 
from C, N, O, and Fe, while Ca4455 is hardly altered by those elements. 
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Fig. 8. — At 4000 K and log g = 2.5, solar-scaled and Ca-enhanced spectra are divided. The index definition of lWorthev et al] JTOOj) for 
the Ca4455 index is depicted with straight lines (solid lines for the index bandpass and dotted lines for the pseudo-continua, respectively). 
Note that the strong Ca line feature is shared both by the index bandpass and the blue continuum that cancel out the calcium effect on 
the Ca4455 index. 
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Fig. 9. — Fe5406 vs. Fe4531 (top) and Fe5015 (bottom). Titanium- and iron-enhanced cases are depicted. Both of these Lick indices 
show a titanium sensitivity equal to that of iron. 
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Fig. 10. — Displacements are shown in diagrams of H/3 as a function of Fe5406 as each chemical element is enhanced by 0.3 dex, except 
carbon, which is enhanced by 0.2 dex. In the lower right panel, solar-scaled chemical mixtures of [Fe/H] = —2.0, —1.0, —0.5, 0.0, and 0.5 
are connected at 1, 2, 5, and 12 Gyr. To guide the eye, the [Fe/H] = 0.5 line is marked with plus signs at given ages. Also in the lower 
right panel, the a-elements enhanced case at fixed [Fe/H] is shown (see text). 
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Fig. 11. — C-, N-, O-enhanced cases are shown from top to bottom in diagrams of H/3 as a function of Fe5406. The same solar-scaled 
chemical mixture grids with a range of metallicity and age are displayed here in the bottom panels as in the bottom panels of Figure 10. 
The left panels show the combination of isochrone and stellar spectral effects, while the right panels display the isochrone effects alone. 
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Fig. 12. — Solar-scaled and Mg-enhanced 12 Gyr integrated spectra near H/3 (4861 A). Note that the Mg-lines make the blue-continuum 
levels significantly lower, which consequently make the H/3 index strengths weaker as we see in the upper right panel of Figure 10. 
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4800 4820 4840 4860 4880 4900 

Wavelength (A) 

Fig. 13. — Solar-scaled and Fe-enhanced 12 Gyr integrated spectra near H/3 (4861 A). Note that Fe-lines make the index bandpass and 
the red-continuum levels significantly lower, which consequently cancel out their effects on the H/3 index strengths as we see in the lower 
right panel of Figure 10. 
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Fig. 14. — Solar-scaled and the a-enhanced 12 Gyr integrated spectra near H/3 (4861 A). Note that the a-enhanced spectra at constant 
Z mostly reflect the Fe depression effect which is the opposite of what is seen in Figure 13 and they affect the H/3 index strengths little as 
we see in the lower right panel of Figure 10. 
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Fig. 17. — Solar-scaled and Fe-enhanced 12 Gyr integrated spectra near H7 (4341 A). The solid lines and the dashed lines are the index 
bandpasses (the broader ones are the A indices, while the narrower ones are the F indices) and the continua regions for the A indices, 
respectively, and the dotted lines are the continua regions for the F indices. Jl-fA B.nd }l'yp have the same blue continuum. Note that the 
iron lines make the red-continuum levels significantly lower, which consequently make the H7 index strengths weaker, especially for the 
broader H7yi index as we see in the upper panels of Figure 15. 
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Fig. 18. — Solar-scaled and Fe-enhanced 12 Gyr integrated spectra near H<5 (4102 A). Symbols are same as in Figure 17. Note that the 
iron lines make the continua levels significantly lower, which consequently make the H(5 index strengths weaker, especially for the broader 
H(5^ index as we see in the upper panels of Figure 16. 
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Fig. 19. — The most conspicuously element-sensitive cases are collectively displayed here. Left: G4300 vs. CN2 (top) and Ca4227 
(bottom). Right: Fe5406 vs. Mg b (top) and H/3 (bottom). 
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Fig. 20. — Similar to Figure 19, but now filled symbols denote isochrone effects alone, while the lines show the combination of isochrone 
and stellar spectral effects. It is clear from this figure that stellar spectral effects dominate the changes in the indices, whereas the isochrone 
effects are comparatively insignificant (see text). 
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Fig. 21. — 2 Gyr-old, scaled-solar spectrum smoothed to 300 kms~^ compared with a 2 Gyr-old, +0.3 dex magnesium-enhanced spectr 
with the same Z in the SAURON spectral range. The spectra are normalized at 4750 A. Note that Mgi (~ 5100 A) and Mg b (~ 5175 
become stronger with magnesium enhancement. 



Stellar Population Models and Individual Element Abundances. II. 



39 



X 

0) 
N 
• I — I 

a 

o 



0.9 



0.8 



0.7 - 



"1 I I r 



"1 r 




solar-scaled 
Fe-enhanced 



J I I L 



4900 5000 5100 5200 
Wavelength (A) 



5300 



Fig. 22. — 2 Gyr-old, scaled-solar spectrum smoothed to 300fcms~^ compared with a 2 Gyr-old, +0.3 dex iron-enhanced spectrum with 
the same Z in the SAURON spectral range. The spectra are normalized at 4750 A. Note that Fe5015, Fe5270, and Fe5335 become stronger 
with iron enhancement. 
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Fig. 23. — 2 Gyr-old, scaled-solar spectrum smoothed to 300 fcms~^ compared with a 2 Gyr-old, +0.3 dex titanium-enhanced spectrum 
with the same Z in the SAURON spectral range. The spectra are normalized at 4750 A. Note that Fe5015 becomes stronger with titanium 
enhancement, while the rest of the indices are mostly intact. 
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Fig. 24. — Solar-scaled, Fe-enhanced, and Ti-enhanced 12 Gyr integrated spectra in the Fe4531 index region. Note that the centroids 
of Fe-enhanced and Ti-enhanced (dotted line) spectra are located to the left and right side of the solar-scaled one, respectively, because of 
different locations of Fe and Ti lines within the Fe4531 index bandpass. 



